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In this paper we have studied a model for self-induced aggregation in Brownian particle incorpo- 
rating the non-Markovian and non-Gaussian character of the associated random noise process. In 
this model the time evolution of each individual is guided by an over-damped Langevin equation of 
motion with a non-local drift resulting from the local unbalance distributions of the other individ- 
uals. Our simulation result shows that colored nose can induce the cluster formation even at large 
noise strength. Another observation is that critical noise strength grows very rapidly with increase 
of noise correlation time for Gaussian noise than non Gaussian one. However, at long time limit the 
cluster number in aggregation process decreases with time following a power law. The exponent in 
the power law increases remarkable for switching from Markovian to non Markovian noise process. 

PACS numbers: 



I. INTRODUCTION 

Formation of large spatial structure and clusters by 
the aggregation of small species joining each other con- 
stitutes a broad area of research in pure and applied 
sciences [T}Q]. General properties of aggregation dynam- 
ics include microphysics of clouds and precipitation [T], 
principle of polymer formation^, different types of eco- 
logical problems^, etc., to mention a few. The most im- 
portant direction of the theory of aggregation dynamics 
leads us to the realm of biological systems in which coop- 
erative activity among individuals usually involves social 
behaviours [H [S]. Still now its underlying mechanism in 
various biological systems remains unknown. A number 
of attempt have been made to established rigorous and 
quantitative basis of the emergence of cooperation among 
individuals [5H5]- The kinship theory is the first and im- 
portant step in this direction which is based on genetic 
arguments^lO . But it does not account correctly the ex- 
planation of cooperation that involves among unrelated 
individuals (individuals having no common genes). On 
the other hand, according to some models of mathemat- 
ical population biology and game theory aggregation is 
the result of short and long range interactions among 
different individuals pTl [T^ . According to these models 
the aggregation dynamics of social biological systems has 
the same statistical basis as that have been used in fluid 
dynamics and condensed matter physics. 

The essential requirements for a general and sim- 
ple theoretical description of aggregation dynamics in 
social communities involve the space of state (where 
each point represents the status of the individuals) and 
strategy IQl [13] (that is the rule according to which the 
individuals players decide to change their status in re- 
sponse to partial and complete information about the 
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action of the other players). Based on game theoretical 
approach, Sigmund and NowakjTl] assumed that cooper- 
ation to work in evolved social systems require the knowl- 
edge of reputation or status of their members (players). 
The reputation is denoted by a dynamical coordinate, S 
known as image score. The image score is assigned to 
each player which is the identity to the other member of 
community and indicates both wether an individual pro- 
vides help and if he/she is worthy of being help. Recently, 
Cecconi et al^^^7\ studied the model suggested by Sig- 
mund and Nowak[T4] to test the hypothesis of emergence 
of cooperation by indirect reciprocity among unrelated 
individuals. In their studv [16,, 17. they assumed nonlin- 
ear Fokker-Planck equation for the population of individ- 
uals with a certain image score. The equation had a non- 
local drift term that characterized the strategy. Based on 
the stochastic differential [151 [TO] equations corresponding 
to a Brownian dynamics with drift induced by the spa- 
tial distribution of other random worker the authors have 
also studied the aggregation dynamics |20|. In this study 
they assumed that the random fluctuations correspond 
to the white Gaussian noise. However, the random fluc- 
tuations in the social communities problem are in general 
non thermal in origin. They may appear as a result of 
complicated inherent dynamics and therefore the noise 
of non thermal origin may be non Gaussian and corre- 
lated in characteristics. The correlated noise may have 
important role in the context of cooperation which is re- 
lated to aggregation dynamics and others since correla- 
tion among the particles increases with increase of noise 
correlation time. Keeping in mind of this aspect we have 
extended the study of the aggregation dynamics of the 
self- induced model [5D] of the social biological system in 
non-Markovian and non-Gaussian limit. Our aim is to 
explore how the dynamics of aggregation to form cluster 
depends on non-Markovian and non-Gaussian properties 
of the noise process. Here it is to be noted that the 
noise of biological (non thermal) origin is non Gaussian 
in character. Recent experimental and theoretical stud- 
ies in neural network and sensory systems [HJ [55] offer 
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strong indication that the noise sources in these systems 
could be non-Gaussian. The noise of biological origin in 
many cases is due to nonlinear dynamics which may be 
correlated and non-Gaussian in character, specifically, in 
the context biological evolution [2T1 [51] . The role of col- 
ored non Gaussian noise in the barrier crossing dynamics, 
the stochastic resonance and complex net work has been 
explored by several authors [24j. 



II. THE MODEL 

Consider a system consists of N individuals which 
change their Xi state according to majority rule. Xi de- 
note the reputation score of z*'* member or the position 
in a possible chemotaxis description or some other am- 
plitude characterizing the role of individual within pop- 
ulation biology framework. We assumed that each indi- 
vidual changes Xi by the following stochastic equation of 
motion. 



xj = v{xj) + 'n{t) 



(1) 



V denotes the drift, which is given by the following ex- 
pression 



v{xi) = A 



w+{x,,t) - W^{Xt,t) 



w+{x,,t) - W^{Xt,t) 
where w± are defined by 



(2) 



iJJ±{x) = ^ Q[±{xj - Xi)] exp {-a\xj - Xi\) (3) 
i 

Q is the unitary step function. The equations ([T]) and ([2| 
reveal that the velocity at which an individual decides to 
move to the left or to the right depends on the difference 
— . The magnitude of w± depends on the expo- 
nential weight with a coefficient a = l/rp, specify the 
range up to which one individual still perceives the pres- 
ence and the influence of the other member of the group. 
So according to the above description aggregation is the 
result from preferably migration of individuals depending 
on the sign of v. The term w+ — in the denomina- 
tor of Eq. ([2]) present the normalization factor and so the 
velocity bounded within [—A, A]. 

We assume that r]{t) in the above equation is colored 
noise which may be of Gaussian or non-Gaussian type 
depending on the situation. The non-Gaussian noise can 
be generated from the solution of the following Langevin 
equation [25J 
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1 d 
T drj 



(4) 



C(t) being a standard Gaussian noise of zero mean and 
its two-time correlation given by 



icmt') = 2S{t - 1') 



(5) 
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FIG. 1: (Color online) This figure refer to system dynamics 
with 100 particles staring from uniformly distributed initial 
conditions in the spacial range [—44,-1-44]. The parameter 
set for sub-figures (a) p = 1, t = 0.01, D = 0.01, (b) p = 
1, r = 0.01, D = 0.5 (c) p = 1, r = 0.01, D = 2.0 (d) 
p = 1, T = 5.0, D = 0.5 (d) p = 1, T = 8.0, D = 2.0 (f) 
p = 1, r = 0.01, D = 0.5. a = 1 and A = 1 for the all 
sub-figures 
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T{p-l) 



ln[l + ai(p-l)r,V2)]. 



(6) 



Here the form for noise 77 allows us to control the de- 
parture from the Gaussian behavior easily by changing a 
single parameter p. D and r are noise parameters related 
to the noise intensity and the correlation time of rj. The 
parameter ai in Eq.(|6]) is defined as 



ai = 



D 



(7) 



Now we consider two different situations. 
Eq.(|4]) becomes 



For p 



V 



V 

T 



at) 



(8) 



i. e., the time evolution equation of Ornstein-Uhlenbeck 
noise process ^27j for which the correlation function 
(77(^)77(0)) decays exponentially 



{r,it)rjm^D/T'eM-t/r) 



(9) 



Thus r is the correlation time of the Ornstein-Uhlenbeck 
noise. In the next step we consider a situation where 
p > 1. For 75 > 1 the stationary properties of the noise 77, 
including the time correlation function, have been stud- 
ied in [26' and here we summarize the main results. The 
stationary probability distribution is given by 



1 



l + ai(p-l)| 



(10) 
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where Zp is the normaUzation factor and given by 



1 + q;i(p — 1 

r(i/(p-i)- 



'I 

-1/2) 



-1 
p-1 



ai(p-l) r(l/(p~l)) 



(11) 



r indicates the Gamma function. This distribution can 
be normalized only for p < 3. Since the above distribu- 
tion function is an even function of 77, the first moment, 
(t/), is always equal to zero, and the second moment is 
given by 



r(5 - 3p) 



(12) 



which is finite only for p < 5/3. Furthermore, for p < 1, 
the distribution has a cut-off and it is only defined for 
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t(1-p) 



(13) 



Finally, the correlation time of non-Gaussian noise r of 
the stationary regime of the process 77 (t) diverges near 
p = 5/3 and it can be approximated over the whole range 
of values of p as 



Tp ~ 2r/(5 - 3p) 



(14) 



Clearly, when — >■ 1, we recover the limit of rji being 
a Gaussian colored noise, the Ornstcin-Uhlenbeck pro- 
cess since in this limit the term in the square bracket of 
Eq.(10l can be written as 



9 9 
1] rj 

l + ai{p- 1)— = exp(ai(p- l)y) 



and therefore Eq.(10l becomes 



with 



^(^7) = ^ exp(-ai77V2) 
^1 



(15) 



(16) 



(17) 



Here we would like to note that Eq.(12) shows that for 



a given external noise strength D and noise correlation 
time r the variance of the non-Gaussian is higher than 
that of the Gaussian noise for p > 1, i. e 

('^p) > (^/'> ■ (18) 
Similarly Eq.(2.14) implies that Tp > r for p > 1. 

III. RESULT AND DISCUSSION 

Based on the above mentioned model we have nu- 
merically investigated the aggregation dynamics. To 
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FIG. 2: (Color online) Number of cluster vs time plot for 
different values of the noise strength to point out the critical 
value of noise strength. The parameters for the plots are 
a = 1, A = 1, p = 4, e = 0.25, r = 1, p = 1 



follow the dynamics of the each individuals present in 
the system we solve N number coupled stochastic equa- 
tions (Eq.(2.1)) along with the equation for noise process 
(2.4) simultaneously using standard Heun's algorithm. N 
refers to the number of individuals present in the system. 
A very small time step(Ai) of 0.01 for numerical integra- 
tion has been used. For the initial coordinates we have 
assumed that at t = all the particles are uniformly 
distributed in the space. However, to exhibit the role 
of noise strength, noise correlation time and other noise 
parameters on the cluster formation dynamics we plot 
coordinates of all the particles vs. time, t in Fig.l. It 
shows that the cluster formation tendency as well as size 
of the cluster increases as the noise strength increases up 
to a critical value (Figs.l(a,b)). After the critical value 
of the noise strength all the individual particles exhibits 
normal Brownian motion (Fig. Ic) . But if we increase the 
noise correlation time then even at high noise strength 
cluster formation is possible. It has been shown in the 
Figs.l(d,e). Thus colored noise can induce the cluster for- 
mation. Further more, Fig.le demonstrates that there 
may be a phase transition from the clustered state to 
the state where particles are uniformly distributed over 
the space if one switches from Gaussian to non Gaussian 
noise(Fig.lf). These observations can be rationalize in 
the following way. The cluster formation is a result of 
cooperation among the particles. By virtue of diffusion 
in the presence of noise the particle which is far apparat 
from the nucleation center of the strong cluster may en- 
ter into the cluster zone. Thus noise accelerates to form 
strong and bigger cluster. But if the noise strength is very 
high then the diffusion dominates over the cooperative ef- 
fect and we observe the phase where particle almost uni- 
formly distributed over the space instead of cluster for- 
mation. Now we come to the point how noise correlation 
time can induce the cluster formation. With increase of 
noise correlation time the variance of the noise decreases 
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and thereby diffusion of the particle suppressed as the 
non Markovian nature of the noise grows. Not only that, 
the noise correlation time strongly effects the drift term 
of the dynamics [571 121] ■ The drift term in the present 
problem accounts the extent of cooperation among the 
particles which leads to cluster formation. Thus it is ap- 
parent in the cluster formation dynamics that the colored 
noise induced nucleation is a result of the extension of 
the cooperation and the correlation among the particles 
as well as suppression of diffusion nature of the particles 
with increase of noise correlation time. Finally, we con- 
sider how the cluster formation is suppressed by the non 
Gaussian noise. For a given noise strength the variance 
of the non gaussian noise is much higher than the Gaus- 
sian noise(see Eq.()). As a result of that for a given noise 
strength the diffusion may dominates over cooperative 
effect for non Gaussian noise and there is no cluster for- 
mation. But for the same noise strength the cluster for- 
mation may be possible for Gaussian noise due to weak 
diffusion compared to non Gaussian noise. Thus there 
may exist phase transition phenomenon if one switches 
from non-Gaussian to Gaussian noise or vice versa. 
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FIG. 4: This plot present the variation of the cluster number 
as a function of the correlation time of the noise for different 
values of the non-Gaussian parameter(p). The other param- 
eter sets are A = 1, a = 1, p = 4, e = 0.25. 
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FIG. 3: (Color online)This figure present the critical pa- 
rameter regime for noise correlation and noise strength for 
different values of the non-Gaussian parameters. Other sys- 
tem parameters for the plots are a — 1, A = 1. 



The above discussion implies that there is a critical noise 
strength(Z?c) above which cluster formation is not pos- 
sible for the given parameter set(Figs. l(a)-l(c)). One 
can determine it numerically from the plot of cluster (it 
is defined as a set of particles within a given cutoff dis- 
tance e. e is called resolution.) number(Afc) 'vs. time (t). 
To demonstrate this we have plotted Nc vs. t in Fig. 2 
for the particle density p = 4.0. It shows that for the 
given parameter set for the Fig. 3 the Dc is very much 
close to 0.70 since cluster number remains almost close 
to its initial value. However, to explore its dependence on 
the noise correlation time we have plotted Dc vs. noise 
correlation time in Fig. 3. It exhibits that the value of 
Dc changes rapidly for the Gaussian noise than the non 
Gaussian one. This is because of the higher variance for 



former than later for the given noise strength D. For 
the same reason the value of D^ in general is higher for 
Gaussian noise than that of non Gaussian noise. It is to 
be noted here that the Fig. 3 is a nice demonstration re- 
garding the dependence of phase transition on the noise 
strength(Z?) and the noise correlation time (t). The pa- 
rameter regime below the dashed curve is corresponding 
to the clustered phase for the Gaussian noise. Similarly 
the points below the solid curve represents this phase for 
the non Gaussian noise. 

In the next step we have investigated the variation of 
cluster number at stationary state with noise correlation 
time(T) and plotted in Fig. 4. It exhibits that the clus- 
ter number rapidly falls with r for the Gaussian noise 
than the non Gaussian noise. As the correlation as well 
as the cooperation among the particles increases with in- 
crease of tau the strong and bigger cluster is formed for 
larger noise correlation time. Because of that the cluster 
number reduces for both the Gaussian and non Gaussian 
noises as tow grows. The slow decrease of cluster num- 
ber for non Gaussian noise compared to Gaussian noise is 
due to higher noise variance (after certain value of noise 
variance it is difficult to form cluster) for former than 
later. 

Finally, we come to consider one of the important as- 
pects of the present study. How the aggregation kinetics 
depends on the noise correlation time? To this aim we 
have investigated the aggregation processes through a set 
of simulation on a system involving 10"* particles. Our ob- 
servation shows (Fig. 5) that both in the Markovian and 
non-Markovian limits cluster number rapidly decreases 
with time at short time regime and after that it slowly 
varies . However, For a detail analysis of the influence 
of the noise correlation and non-Gaussian parameter on 
the rate of aggregation processes we set an approximate 
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FIG. 5: The plot presents the variation of the cluster number 
as a function of time at long time limit for different values of 
the noise correlation time (r) and the non-Gaussian parame- 
ter(p). The other parameter sets are A = l, a = 1, p =1 e = 
0.25. In the inset same plot which covers both both the short 
time as well as long time limits. 

power law of the following form. 

N,{t)^t-' . (19) 

This type of algebraic decay law is valid for all t and 
p values provided t is long enough (see the inset of the 
Fig.5). Prom Fig. 5 it is surprising to note that the expo- 
nent z in the power law increases more than fifty percent 
compared to Markovian case for r = 1.0. Thus it puts 
further evidence to consider that the noise correlation 
time enhances the cooperative effect in the aggregation 
dynamics through modification of drift term as well as 
reducing the noise variance for the large noise strength 
case. Before going to leave this section we mention that 
the exponent z is little bit smaller for colored non Gaus- 
sian noise compared to Gaussian colored noise as a result 



of higher noise variance for former than later for the given 
noise strength. 



IV. CONCLUSION 

Based on the numerical simulation of stochastic dy- 
namics associated with colored non-Gaussian noise we 
have studied the aggregation kinetics of an self-induced 
model. Such type of self-induced model is used to deals 
with a certain class of problems in social biological prob- 
lem. In this modelled, each individuals represented by a 
Brownian particle, which undergoes a drift velocity de- 
pending upon the population unbalance perceived by a 
single individual between left and right. Through a set of 
numerical simulation we have examined the effect of the 
noise correlation and non-Gaussian character of the noise 
in the self-induced aggregation dynamics. Our main ob- 
servation includes the following points: 

(i) The colored noise can induce aggregation even at 
large noise strength (Fig.l). The aggregation may dis- 
appear if one switches from Gaussian to non Gaussian 
noise. 

(ii) The critical noise strength rapidly increases with 
noise correlation for non Gaussian noise than the non 
gaussian one. 

(iii) The cluster formation kinetics follows a power law 
for the variation of cluster number with time like, N^it) ~ 

at long time. The exponent z remarkable increases 
for non Markovian case. 

We hope that our observations will be useful to under- 
stand the aggregation processes in the social biological 
system. 
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